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The masses and decay constants of the ground statě heavy XQ2 (Q = b, c) tensor 
. mesons are calculated in the framework of the QCD sum rules approach. The obtained 

results on the masses are in good consistency with the experimental values. Our 
■ predictions on the decay constants can be verified in the future experiments. 

(N ; PACS number(s): 11.55.Hx, 14.40.Pq 

(N ' 
O 
O 



X 



*e-mail: taliev@metu.edu.tr 
^permanent address:Institutc 
íe-mail: kazizi@dogus.edu.tr 
§e-mail: savci@metu.edu.tr 



of Physics, Baku,Azerbaijan 



1 Introduction 



During last few years very exiting experimental results are obtained in the charm and 
beauty meson and baryon spectroscopies [1]. Recent CLEO measurements on the two 
photon decay rates of the even-parity, scalar ++ , Xř>( c )o an d tensor 2 ++ ,Xb(c)2 states ([1,2] 
and references therein) were motivation to investigate the properties of these mesons and 
their radiative decays. 

In the present work, we calculate the mass and decay constants of the ground statě heavy 
bottomonium, Xb2^-P) and charmonium, xaO-P) tensor mesons with I G (J PC ) = + (2 ++ ) 
in the framework of the QCD sum rules approach. QCD sum rules approach as a non- 
perturbative approach is one of the most powerful and applicable tools to hadron physics 
and can play an important role in calculation of the characteristic parameters of the hadrons 
(for details about this method and some applications see [3,4]). Note that the mass and 
decay constant of the strange tensor i^|(1430) with quantum numbers I(J P ) = 1/2(2+) 
have been calculated in [5] in the same framework. These parameters for light unflavored 
tensor mesons have also been calculated in [6]. The obtained results for the masses and 
decay constants are used in calculation of the magnetic dipóle moments of the light tensor 
mesons using the QCD sum rules method in [7]. 

The páper is organized as follows: in next section, sum rules for the mass and decay 
constant of the ground statě heavy quarkonia, XQ2 tensor mesons are derived in the context 
of the QCD sum rules method. Section 3 is devoted to the numerical analysis of the mass 
and decay constants as well as the comparison of the obtained results on the mass with the 
experimental values. 

2 Theoretical Framework 

In this section, we obtain the sum rules for the mass and decay constant of the heavy 
XQ2O-P) tensor meson in the framework of the QCD sum rules approach. For this aim we 
consider the following correlation function 



where, is the interpolating current of the xq 2 (1P) tensor meson and T is the time 
ordering operátor. The explicit form of the current j M „ creating the ground statě heavy 
tensor %q 2 (1P) statě with quantum numbers I G (J PC ) = + (2 ++ ) from the vacuum can be 
written in the following form: 



where Q stands for heavy b or c quark and the (x) represents the derivative with respect 
to four-x acting on left and right, simultaneously. This two-side covariant derivative is 
defined as: 




(1) 



\Q(x)^ V v (x)Q(x) + Q{x)^ v V» (x)Q(x) 



(2) 




(3) 
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where, 

^»(x) = -Ě,(x)-i 9 -\ a Al(x), 

^^x) = ^^x)+i^A%{x). (4) 

In the above relations, the A a are the Gell-Mann matrices and A a {x) are the external 
(vacuum) gluon fields , which can be expressed directly in terms of the gluon field strength 
tensor using the Fock-Schwinger gauge, x^A a (x) = 0, as the following way: 

A;(x) = £ daax^ax) = ^G^(O) + ^x v x p V v G^(0) + ... (5) 

Since the current contains derivatives with respect to the space-time so we will consider the 
two currents at points x and y. After doing calculations and applying the derivatives, we 
will set y — 0. 

It is well known that in the QCD sum rules approach, the correlation function in Eq. 
(1) is calculated in two different ways. The physical or phenomenological part, which is 
obtained in terms of the hadronic parameters such as mass and decay constant inserting a 
complete set of the states owing the same quantum numbers as the interpolating current 
The theoretical or QCD part, which is calculated in terms of the QCD parameters 
such as quark masses, quark condensates and quark-gluon coupling constants etc. The 
correlation function in this part is calculated in deep Euclidean region, q 2 <C 0, via operátor 
product expansion (OPE). The short distance effects are calculated via the perturbation 
theory, whereas the long distance contributions, which are non-perturbative effects are 
parameterized in terms of the quark-quark, gluon-gluon and quark-gluon condensates. The 
sum rules for the observables (masses and decay constants) of the ground statě x<22(l-P) 
meson are obtained equating both representations of the correlation function, isolating the 
ground statě and applying Borel transformation to suppress the contribution of the higher 
states and continuum through the dispersion relation. 

To proceed first we calculate the phenomenological part. Inserting a complete set of 
intermediate statě, XQ2^-P) to time ordering product in Eq. (1), and performing integrál 
over x we get: 

TT (0 I 3j»M I XQ2)(XQ2 I j a p(0) | 0) ( . 

ííuu,a/3 — 5 5 T • • • , [V) 

m XQ2 ~ 1 

where • • • denotes the contribution of the higher states and continuum. From the above 
relation, it is clear that we need to know the matrix element, (0 | j^(0) | XQ2), which can 
be parameterized in terms of the decay constant, f XQ2 : 

(0 | J^(0) | XQ2> = f X Q2 m XQ ^^, (7) 

where e^ u is the polarization tensor of XQ2 meson. Using Eq. (7) in Eq. (6), we obtain the 
following finál representation of the correlation function in phenomenological side: 

fl 02 rrň fi 1 
= 2 _ 2 1 9(sWS^/3 + flWJflw) \ + other structures + (8) 

m XQ2 1 l Z ) 
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where, the only structure which contains a contribution of the tensor meson has been kept. 
In calculations, we have performed summation over the polarization tensor using 

1„ „ 1„ „ 1, 

3 

where, 



= + M^. (10) 

The next step is to calculate the theoretical or QCD side of the correlation function 
in deep Euclidean region, g 2 < . Using the explicit expression for the tensor current 
presented in Eq. (2) inside the correlation function shown in Eq. (1) and contracting out 
all quark pairs using the Wick's theorem, we obtain the following expression for the QCD 
side: 

IW = \ J d A xé^ |Tr [s Q (y - x)^ V v (x) Vp (y)S Q (x - j/) 7a ] 
+ [floa] + [fO/í] + [j9«a,!/e/i]J. 



(11) 



To proceed we need to know the heavy quark propagátor, Sq(x — y). This propagátor has 
been calculated in [81: 



Sq{x-y) = S[ ree (x -y) - ig s J 



d A k 



-ik(x-y) 



dv 



fé + rriQ 
(m 2 Q — k 2 ) 2 



G^[v(x-y)]a 



fj,i/ 



+ 



1717 



Mx, - Vll )cr lv 



;i2) 



where, 



and K n {z) being the modified Bessel function of the second kind. The next step is to use 
the expression of the heavy propagators and perform the derivatives with respect to x and 
y in Eq. (11). After setting y — 0, the following finál expression for the QCD side of the 
correlation function in coordinate space is obtained: 

(^) 4 / d ^ qX {[ r M^/3] + \P «] + [" tA + \P a, v O »]} , (14) 



n, 



where, 

r 



% 
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- IrriQgavgpiJCíK-i + 2\m Q x a x l/ g^ - g av g^ + g afl gp v + g a pg^jíC 2 

- 2m 2 Q x Q x v g fiil K x K z - 2m Q g au (2x p x^ - x 2 g f3 ^j /C 2 /C 3 

+ 2m 2 Q x a x u (2x p x^ - x 2 g Pl }j K\ - 2m 2 Q x a x v {^xpx^ - x 2 g P! ^ £ 2 /C 4 } 

+ nonperturbative contributions , (15) 
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and 

K n = *"(™g^) . (16) 

In the present work, we calculate the contributions of the heavy quark and gluon con- 
densates in nonperturbative part of the correlation function in QCD side. After a simple 
calculation we obtain for the heavy quark condensate (for the coefficient of the aforemen- 
tioned structure) 



rrin 



(QQ) 



2tf-m 2 Q ) 

Using the well-known relation between the heavy quark and the gluon condensates 

these two nonperturbative contributions can be written in terms of gluon condensate con- 
tribution. Numerical analysis shows that, taking into account quark condensates decreases 
gluon condensate contribution about 15%. 

Few words about the neglected dimension two operátor in the operátor product expan- 
sion are in order. The term proportional to l/g 2 introduced in [9] is a phenomenological 
parametrization of the higher order contributions to the perturbative series. In other words, 
this term can appear when considering any types of correlation functions where the per- 
turbative series are not zero. Obviously, this term vanishes when considering the difference 
of the correlators induced by vector and axial vector currents, VV-AA in the chiral limit, 
m q = (for more details see [10]). In the present work, we neglect this term because we 
work to leading order in a s . 

Now, we apply the Fourier transformation to the QCD side of the correlation function 
to get its expression in momentům space. The next step is to select the structure which 
gives contribution to the tensor statě from both sides of the correlation function, equate 
the coefficient of the selected structure from both sides and apply the Borel transformation 
to suppress the contribution of the higher states and continuum. After lengthy calcula- 
tions, finally we obtain the following sum rules for the decay constant of the heavy tensor 
quarkonia: 

- 2[4M 2 - s - s(u)]u + [2m 2 Q + 4M 2 - s - s(u)]u 2 } + I(M 2 ) (—G 2 ^ , 

(17) 



where, 



u + 



1 



(18) 
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and the explicit expression of the function I(M 2 ) is quite lengthy and therefore we do not 
present it. 

In the above sum rules, M 2 is the Borel mass parameter, s is the continuum threshold 
and N c = 3 is the color factor. The mass of the heavy tensor meson is also obtained 
applying derivative with respect to — to the both sides of the sum rules for the decay 
constant and dividing by itself, i.e., 



( r° , r , e ~ s/ v - s s(u)\ f n 2 3 r ~ 

- 2[4M 2 - s - s(u)]u + [2m 2 Q + 4M 2 - s - s(u)]u 2 } 

+ 1; * r * ^^^{^ -)■} - (v- 2 ) ) 



2[4M 2 - s - s(u)]u + [2m£ + 4M 2 - s - s(u)]u 2 } + I(M 2 ) (^—G : 



(19) 



3 Numerical analysis 

In this section, we numerically analýze the sum rules for the mass and decay constant 
of the ground statě tensor quarkonia. Some input parameters entering the sum rules are 
quark masses, m b = (4.8 ± 0.1) GeV, m c = (1.46 ± 0.05) GeV [4] and gluon condensate, 
(0|-a: s Gr 2 |0) = (0.012±0.004) GeV A . It should be noted that recent analysis of experimental 
data leads to the (0\±a s G 2 \0) = (0.005 ± 0.004) GeV 4 for the gluon condensate [11]. For 
conservative estimation in numerical analysis, we also take into account the value of gluon 
condensate (0\±a s G 2 \0) = (2.16 ± 0.38) x 10~ 2 GeV 4 which follows from sum rules of 
e + e~ — > I — 1 hadrons [12] and heavy quarkonia [13-15]. Few words about quark mass are 
in order. The aforementioned masses are the pole masses for the quarks. Using the four 
loop results for the vacuum polarization operátor in [16], the running masses of the charm 
and beauty quarks, m c (3GeV) = (0.986 ± 0.013) GeV and m b (m b ) = (4.163 ± 0.016) GeV 
are obtained. These improved values of the running masses of charm and beauty quarks 
as well as wide range of gluon condensate are used in numerical calculations. To obtain 
more reliable results for the mass and decay constant of the heavy tensor meson, we will 
also take into account a more realistic error coming from the range spanned by the pole 
and running quark masses as well as the range for the value of the gluon condensate. 

From the sum rules for the decay constant and mass it is clear that they contain also 
two auxiliary parameters, continuum threshold s and Borel mass parameter M 2 . The 
standard criteria in QCD sum rules is that the physical quantities should be independent 
of these mathematical objects, so we should look for working regions for these parameters 
at which the masses and decay constants practically remain unchanged. To determine the 
working region for the Borel mass parameter the proceduře is as follows: the lower limit 
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of M 2 is obtained requiring that the higher states and continuum contributions constitute, 
say, 30% of the total dispersion integrál. The upper limit of M 2 is chosen demanding that 
the sum rules for the decay constants and masses should be convergent, i.e., contribution 
of the operators with higher dimensions is small. As a result, we choose the regions: 
8 GeV 2 < M 2 < 20 GeV 2 and 4 GeV 2 < M 2 < 7 GeV 2 for the Borel mass parameter. 
The continuum threshold Sq is not completely arbitrary but it is correlated to the energy 
of the first exited statě with quantum numbers of the interpolating current. Our numerical 
results are in consistency with this point and show that in the interval ijn XQ2 + 0.4) 2 < 
Sq Q2 < {m XQ2 +0.7) 2 , the results are practically insensitive to the variation of this parameter. 
Here we would like to make the following remark. It is shown in [17] that the continuum 
threshold Sq can depend on the Borel mass parameter. Therefore, the standard criteria, 
námely, weak dependence of the results on variation of the auxiliary parameters does not 
provide us realistic errors, and in fact the actual error should be large. Following [17], in 
the present work we will add also the systematic errors to the numerical values. 

Our numerical analysis on the masses and decay constants leads to the results presented 
in Table 1. The quoted errors in our predictions are due to the variations in the continuum 





Present Work 


Experiment [1] 




(9.90 ±2.48) GeV 


(9.91221 ±0.00057) GeV 


m Xc2 


(3.47 ±0.95) GeV 


(3.55620 ± 0.00009) GeV 


fxb2 


0.0122 ±0.0072 




fxc2 


0.0111 ±0.0062 





Table 1: Values for the masses and decay constants of the tensor mesons xq 2 - 

threshold and Borel parameter, uncertainties in quark masses and wide range of the gluon 
condensates presented at the beginning of this section as well as the systematic errors. 
The results presented in Table 1 show a good consistency between our predictions and the 
experimental values [1] on the masses of the ground statě heavy, Xb2 and, Xc2 tensor mesons. 
Our predictions on the decay constants can be verified in the future experiments. 
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